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In  a  number  of  areas  of  science  and  technology  it  is  necessary  to 
deal  with  rhythmic  processes,  such  as  the  parasitic  modulation  of  shot 
fluctuations  in  radio-electronic  equipment  [8,9],  changes  in  the  levels 
of  stability  of  oscillatory  parameters  of  biological  systems  [12], 
modulation  of  an  acoustic  field  in  a  layer  of  water  due  to  heaving  of 
its  surface  [7,  13],  change  in  meteorological  parameters  under  the 

action  of  tidal  forces  [14],  etc.  Many  investigators  [1,  11,  12]  note 

the  advisability  of  using  methods  of  the  theory  of  stationary  random 
processes  for  analyzing  such  essentially  nonstationary  phenomena  as 
biological,  etc.  rhythms,  and  stress  the  necessity  of  developing 
statistical  methods,  based  directly  on  the  properties  of  this  type  of 
nonstationary  processes  and  having  an  explicit  interpretation.  The 
frequently  used  additive  model  -  the  determinate  periodic  function 
^(t)  on  a  background  of  stationary  noise  V{t)-  j [t)+  V(t) ,  analyzable 
by  methods  of  the  theory  of  processes  with  stationary  pulsations,  does 
not  describe  the  main  property  of  rhythmic  phenomena  -  change  in  the 
parameters  of  the  process,  analogous  to  pulsations  (of  modulation).  A 
multiplicative  model  of  the  type  of  modulation  -  ■^(t),  where 

(t)  -  stationary  process,  although  it  describes  this  feature  it  is 
not  general.  Therefore  the  need  arose  for  developing  a  general 
probability  model  of  rhythmic  phenomena  and  on  this  base  to  develop 
effective  methods  for  analyzing  them.  Although  in  the  literature  [11] 
the  important  role  of  a  study  of  periodic  nonstationary  random 
processes  is  noted  in  connection  with  an  analysis  of  rhythmic 
phenomena,  this  thesis  is  in  no  way  substantiated,  and  the  theory  of 
this  class  of  processes  is  found  in  the  initial  phase  of  development, 
to  say  nothing  of  its  applications.  In  work  [8]  only  the  process  of  a 
type  of  modulation  is  studied.  Taking  this  into  consideration,  the 
author  together  with  K.S.  Voychishin  made  an  attempt  at  developing  and 
substantiating  a  general  model  of  rhythmic  phenomena.  The  idea  of 
creating  such  a  model  as  a  base  for  developing  methods  of  analyzing 
them,  the  use  for  this  purpose  of  periodically  correlated  (PK)  random 
processes  (SP)*,  and  the  proof  of  the  necessary  theorems  belong  to  the 
author,  and  for  substantiation  of  the  model  the  results  from  processing 
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the  materials  of  observations  c£  the  bioactivity  of  loaches  and 
meteorological  elements  which  were  obtained  by  K.S.  Voychishin  were 
used  considerably. 

*  Footnote.  Here,  just  as  in  work  [5],  we  will  use  the  term 
"periodically  correlated"  in  the  sense  of  "periodically  nonstationary 
of  the  second  order"  (or  in  a  wide  sense,  if  an  analogy  with  stationary 
processes  is  used).  This  class  of  processes  is  sufficient  for  the 
theory  being  developed.  Periodic  nonstationary  in  a  narrow  sense  is 
defined,  as  usual,  through  the  periodic  nature  of  n-di mensional 
distributions  of  probabilities  [11].  [End  of  footnote] 

For  creating  a  meaningful  theory,  in  contrast  to  mathematical,  it 
is  not  sufficient  just  to  select  axioms,  it  is  also  necessary  that  they 
represent  significant  sides  of  the  phenomenon  for  which  the  theory  is 
being  constructed.  Therefore  in  work  [2]  on  the  basis  of  an  analysis 
of  tracings  of  rhythmic  phenomena  it  is  postulated  that  if  a  phenomenon 
takes  place  in  an  unchanged  mode,  then  it  is  reasonable  to  accept  the 
following  as  its  significant  characteristics:  a)  noise-like  nature, 
i.e.,  the  fundamental  impossibility  of  their  representation  in  the  form 
of  a  complex  determinate  function;  b)  harmonizani 1 i ty  -  separability 
into  simple  harmonic  components;  c)  nonstationary  state  of  the  type  of 
modulation,  characterized  by  a  periodic  change  in  dispersion,  and  this 
means  in  covariation. 

From  these  postulates  it  follows  [6]  that  in  the  correlation 
theory  a  periodically  correlated  random  process  (PKSP)  should  be  a 
general  model  of  a  rhythmic  phenomenon. 

The  purpose  of  this  article,  just  like  that  of  article  [6], 
developing  the  ideas  of  work  [5]  and  supplementing  work  [2],  is  a 
further  study  of  the  structural  properties  of  PKSP  as  a  model  of 
rhythms,  but  in  contrast  to  article  [6]  the  mathematical  expectation  of 
PKSP,  generally  speaking,  is  considered  different  from  zero. 

Def inition  1 .  We  will  call  the  random  process  ( SP )  centered  if 
its  mathematical  expectation  is  equal  to  zero. 

Then  from  a  presentation  of  PKSP  through  stationary  components 
(SK)  [6] 


and  the  uniqueness  of  the  expansion  into  a  Fourier  series  follows  the 
validity  of  the  theorem. 


2 


Theorem  1.  A  PKSP  is  centered  when,  and  only  when,  its  stationary 
components  are  centered. 

Theorem  2.  If  the  SK  PKSP  is  uncorrelated  and  centered,  with  the 
exception  maybe  o£  zero,  then  it  degenerates  into  into  a  stationary  SP, 
if  the  SK  are  only  uncorrelated,  then  into  SP  with  stationary 
pulsations  and  a  periodic  mean. 

Actually,  from  correlation  [5] 

«  »»  “ 

b^{t,u)=  V  B^{u)e  =  V  ^  [?} 

Oft  k^-’-Ca  Iwm — OO 

(where  r j  {  •  )  -  mutual  covariation  of  k-  and  j-th  SK  under  the 

condition  of  noncorrelation  of  SK  r  j  {  u  )  =  r  j,;  {  u  )  S  j,k'  from  which 
^1  ,  1-k  (  u  )  =ri  (  u  )  5  j,;  _  0  )  we  obtain 

^  ^  V  r,  (u) « ' ^  (u). 

where  Bg ( • )  -  covariation  of  a  stationary  approximation  to  PKSP  [5], 
i.e.,  SP  ^  (t)  is  stationary  relative  to  covariation.  And  for  the 
mathematical  expectation  from  formula  (1)  we  find 

• 

m|(0=  £  «»«  . 

‘'km—m 


where  mj^  =  E  -  in  a  general  case  (for  all  mj^^O)  a  periodic 

function  with  a  period  of  correlation. 

Theorem  3.  The  non-correlatedness  of  odd  SK  PKSP  with  even  and 
paired  together  is  sufficient  so  that  its  covariation  would  possess 
half  a  period. 

Actually,  under  the  conditions  of  the  theorem 

jo,  if  (  j  V 
1,  if  {  j /\  k  )  are  even, 

•»  SB 

and  then  in  equality  (2)  ^*(“)=  £  identically 
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equal  to  zero  in  the  case  of  odd  k. 


Cor o 1 1  ary  1 .  Under  the  conditions  of  the  theorem  a  FKSP  possesses 
a  half  a  period  if  its  odd  SK  are  centered. 

Corollary  2 .  The  absence  of  SK  is  sufficient  so  that  the  FKSP 
would  possess  half  a  period. 

Corollary  3 .  Centered  SK,  non-correlated  pairwise  and  with  all 
the  remainder,  can  be  related  to  zero. 

When  all  the  SK  are  linear  transforms  of  one  of  them  {for  example, 
this  will  happen  when  the  PKSP  is  obtained  by  conversion  of  a 
stationary  SP  by  a  system  with  periodically  changing  parameters): 

5*  (0  =  A,l  (0  -  J  (X)  dZ(\). 


where  •)  -  frequency  characteristic  of  the  corresponding  invariant 

filter  ,  then 


5(0  = 


and  the  correlation  components  (KK) 


S  t(M-/ 


and  the  covariations  of  components 

r*.(u)= 


From  here  it  is  evident  that  the  components  will  be  non-correlated  if 

In  the  particular  case  under  consideration  it  is  directly  evident 
that  a  situation  is  possible  when  B|^(0)=0  and  Bi^(u)BjfcrO  in  the  case  of 
usfsG  •  For  this  it  is  sufficient  to  assume  that  the  j-th  SK  is  the  k-th 
Hilbert  transform  ,  i.e.,  TP*  {X  and  Ipj  (X  )=-i  sgn  %  '  since 

^./(w)  “2  J$inXdSj(X). 


Example.  Assume 


’1(0  =  5(0  »'■''  + 5(0*"*'^. 

where  V  -  denotes  the  conversion  to  the  Hilbert  transform,  then 
taking  mtc  account  we  obtain 

(^>  “)  =  2rj  («)  cos  A  u  +  (u)  ^  where  ^,  (“)  = '^  (— u)  -  from  here 

B2  =  0.  The  dispersion  of  this  process  o* (0  =  2^^ (0)  =  2oj  =  .  If  the 

process  ^  (t)  is  noncentered  and  E^(t)=m,  then  my^(t)=2m  cos  /At. 

In  the  case  when  the  converter  is  a  diffeiential  operator  with 

M 

periodic  coefficients  {a'  (0. 1  =  i.N)  ,  then  ♦*  (M  =  2]  '  where  aj 

M 

the  k-th  coefficient  of  the  Fourier  function  al(*).  when  only  the 
actual  process  in  the  steady-state  mode  are  considered,  then  odd  1  are 
absent  and  the  period  of  cor  relatedness  of  the  process  will  be  equal  to 
half  of  the  period  of  change  of  the  parameters  of  the  converter.  Then, 
having  designated  the  product  of  the  series 

'■ww-  £  «;(■!)'  2  i;  [  i  . 


for  the  covariation  SK  we  find 

r,/(«)=  Je'“^P*,(X)dS(X). 


If  we  designate  j+ (j.)  =  i ^ an  even  part  of  spectral  density  and 
take  into  account  the  evenness  of  functions  a}  {■)  ,  then 

'u.v  (0)  =  2  J  (X*) ,+  (xj  jx. 

0 


since  the  integral  with  an  odd  part  is  equal  to  zero.  From  here  the 

o» 

initial  value  of  KK  according  to  formula  (2)  (0)  =  J  (X*)  j"*"  (X)  ifX  /  where 

Q«(X’)»2V  ,.(XS)  For  example,  for  operators  of  the  second  and 

\  *  / 

t<m0 
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fourth  orders 


«  «o 


It  follows  from  this  that  the  values  of  KK  in  zero  only  in 
exceptional  cases  will  be  zeros,  i.e.,  for  steady-state  processes, 
representing  transforms  of  stationary  processes  with  the  help  of 
differential  operators  with  periodic  coefficients,  as  a  rule,  the 
dispersion  and  mathematical  expectation  (if  the  processes  are  not 
centered)  possess  a  precise  period,  equal  to  the  period  of 
CO r r e 1  a t edness  and  are  half  of  the  period  of  change  of  the  operator 
coefficients.  Thus  the  PKSP,  the  result  of  a  stationary  disturbance  of 
a  system  with  periodically  changing  parameters,  will  possess  a  period 
which  is  half  the  period  of  change  of  its  parameters. 

The  null  covariation  component  describes  the  distribution  of 
energy  by  harmonics  which  make  up  the  PKSP  [6],  and  therefore  its  role 
in  the  characteristics  of  the  process  is  not  at  all  different  from  the 
role  of  other  KK,  which  is  also  evident  from  the  following  theorems. 

Theorem  4.  In  order  that  the  zero  KK  be  identically  equal  to  zero 
it  is  necessary  and  sufficient  that  it  be  equal  to  zero  in  zero. 

This  follows  directly  from  the  fact  that  Bo(')  is  a  covariation  of 
a  stationary  approximation  to  PKSP,  and  the  covariation  of  the 
stationary  process  satisfies  the  condition  r »  ( 
this  case  gives  the  inequality  /Bq  ( u  )/4:Bo  (  u  )  . 

From  formula  (2)  we  find  that  fl,(0)  =«  r/(0)  ,  and  since  all 

r,ld)  =  o*>0  as  dispersions  of  stationary  processes,  then  the  equality 

Bo(0)=0  is  possible  only  in  the  case  when  all  0^  =0.  But  since  from 

«• 

this  same  formula  (2)  and  the  mutual  covariation  of 

stationary  processes  satisfies  the  inequality  \ftk{u)\*<ri{0)r^{0)  -  then 

from  here  follows  the  conclusion,  which,  in  view  of  its  importance,  we 
will  formulate  as  a  lemma. 

Lemma  1.  If  Bo(0)=0,  then  all  KK  are  identically  equal  to  zero. 


u).<r  -  (0)=  Q-?  ,  which  in 
5  S 
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Actually,  since  when  Bq{0)=0  pulsations  are  absent,  the  process 
degenerates  into  a  determinate  function.  Thus  the  following  theoiem 
takes  place . 

Theorem  5.  Equality  to  zero  of  the  value  in  zero  of  a  zero  KK  is 
necessary  and  sufficient  so  that  the  PKSP  would  degenerate  into  a 
determinate  function  (equal  to  its  average). 

The  sufficiency  of  the  condition  follows  from  formula  (8),  lemna  1 
and  the  definition  of  dispersioi  .  And  the  necessity  is 

evident  from  the  fact  that  the  dispersion  of  pulsations  of  a 
determinate  function  is  equal  to  zero.  Then  from  the  property  of 
covariation  I  "I"  “)  (0  follows  that  |6^(f,u)|<;0  ,  and 

since  ^*  (“)  1  ^  f  j  I  I  then  Bj^(u)=0  in  the  case  of  all  k  and  u. 

0 

Def inition  2 .  As  normalization  of  SP  we  will  name,  in  an  analogy 
with  normalization  of  random  variables,  the  transition  from  process 
t  ( t )  to  process 


l(t)  S(0-«i(0  (3) 

Its  covariation 


1(0 


(/.  U)  = 


6j  (1.  u) 

Oj(X  +  U)Oj(<)’ 


(4) 


and  the  average  and  dispersion,  naturally,  are  defined  as  (0  =  Ot  (f)  =  I 
in  the  case  of  all  t.  Therefore  the  unique  characteristic  of  rhythm  of 
a  normalized  PKSP  can  only  be  covariation.  Since  formulas  (3)  and  (4) 
have  meaning,  when  (Tj  ( t )  *ft*0  for  all  t,  then  from  formula  (2)  it  follows 

- 

t^at  the  function  (^  +  «)  (0  =  V  <  V  B^{0)  e  ^ 

is  periodic  with  respect  to  t  with  a  period  T  and  is  non - van i sh i ng . 


Then  based  on  the  Wiener  theorem  [3]  the  function  (o^  (/ 4*  Oj  (01  * 

is  expanded  into  an  absolutely  convergent  trigonometric  series  of  the 
same  order . 
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Since  the  product  of  the  periodic  functions  of  one  period  is  a  periodic 
function  of  the  same  period,  then  this  proves  the  validity  of  the 
theorem. 


Theorem  6.  The  periodic  correlatedness  of  an  SP  is  invariant  with 
respect  to  its  normalization. 


In  particular, 
constant  dispersion 


6-- (/,«)  = 


5n 


the  covariation  of  a  normalized  PKSP,  possessing  a 
(evidently  (t)=Bo(0)),  will  be 

It  is  advisable  that  normalization  of  processes 


by  applied  in  the  case  of  a  comparative  study  of  the  rhythm  of 
diverse-scale  physical  and  other  parameters. 


Assume  now  that  (t)  -  a  stationary  SP  with  the  covariation 

«• 

R^(-),  and  fit)  »  ^e‘*^dF(k)  ,  where  F(*)  -  certain  complex  measure  of 

limited  variation,  then  the  process 


S(0  =  Ti(0/(0 


(5) 


is  harmonizable  [10]  and  its  function  of  covariation 


6,  (/.«)=/?,(«) /(/  +  «)/(/). 


(6) 


and  dispersion 


{'/ J 


Then  the  normalized  function  of  covariation 

where 
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Dejinition  3^  A  process  of  the  type  (5)  we  will  call  a  process  of 
the  type  of  simple  modulation. 

If  the  function  "f  i  '  )  is  periodic  and  {-j.-}  -  its  Fourier 

coefficients,  then  the  KK  of  the  process  ^  ( t )  will  be 

The  following  theorems  stem  from  the  previous  formulas. 

Theorem  7.  A  process  of  the  type  of  simple  modulation  with  an 
actual  periodic  modulating  function  is  a  PKSP  which  has  been  normalized 
t.'.'  stationary. 


Actually  under  the  conditions  of  the  theorem  b  t  ,  u  )  =  u  )  . 

Theorem  8.  A  symmetry  of  the  third  kind  of  a  modulating  function 
of  an  SP  of  the  type  of  simple  modulation  is  a  necessary  and  sufficient 
condition  of  the  fact  that  the  period  of  its  dispersion  would  be  half 
the  period  of  the  correlated  state. 


The  condition  of  symmetry  of  the  third  kind  is  fulfillment  of  the 

equality  /  ^  ^  )  = — /(Oi*^  case  of  all  t,  therefore  the  function 

7 

//(•)/  possesses  a  period  and  confirmation  of  the  theorem  follows 

from  formula  (7).  ^ 


Theorem  9.  Simple  modulation  is  necessary  and  sufficient  for  the 
conversion  of  white  noise  into  nonstationary  white  noise,  and  if  the 
modulating  function  is  periodic,  then  this  is  PK  noise  [6]. 

The  necessity  follows  from  the  definition  of  nonstationary  white 
noise  and  the  fact  that  the  negative  function  is  a  square  of  the 
modulus  of  the  function  g { t ) = / f ( t ) ^ 2^0 ,  and  the  sufficiency  is  found 
from  formula  (6)  and  the  property  of  the  -function: 

5(x-a)  /^  (  X  )  =  5  {  x-a  )  f  {  a  )  . 

From  the  analysis  made  of  the  PKSP  it  is  evident  that  the  precise 
period  of  dispersion  depends  on  the  values  in  the  zero  of  the  KK,  and 
the  precise  period  of  the  average  -  on  the  centrality  of  the  SK .  Since 
in  a  general  case  these  properties  are  independent,  then  the  period  of 
dispersion  and  the  period  of  the  average  can  be  used  for  estimating  the 
precise  period  of  correlatedness  of  the  PKSP.  In  this  case  the 
dispersion  is  more  significant  than  the  average,  which  frequently  is 
considered  equal  to  zero  [4].  The  most  significant  characteristic  is 
the  function  of  covariation,  describing  the  correlatedness  of  the 
harmonic  components  of  the  PKSP,  and  its  period  is  determined  by  the 
KK  . 
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In  particular,  for  a  normalized  PKSP  the  function  of  covariation  is  the 
unique  characteristic  of  rhythm. 

It  is  evident  from  formula  (1)  that  when  all  SK  are  proportional 
to  one  of  them,  which  can  be  considered  as  zero,  i.e.,  when 
5k(t)=o<(^lo('t),  where  (otk)  -  complex  numbers,  then  the  PKSP  turns 
into  a  process  of  the  type  of  simple  modulation 

5(0  =  5o(0  ^  '  =  ^(0/(0. 


If  all  the  SK  are  non-cor related  and  non-centered ,  then  the  PKSP  turns 
into  an  additive  model  -  a  periodic  determinate  function  on  a 
background  of  stationary  noise: 

5(0=  S  m*e“~'  +  n(0. 


where  (t)  -  stationary  SP  with  a  null  average  and  the  covariation 

=  S  r^(u)e  ^ 

It  is  evident  from  the  considerations  presented  that  the  idea  of 
isolation  of  the  determinate  component  of  the  rhythmic  phenomenon  and 
the  study  of  the  stationary  residual  turns  out  to  be  inconsistent, 
since  only  in  an  additive  model  can  it  have  direct  meaning.  In  this 
case  the  stationary  residual  cannot  describe  the  rhythmic  nature,  and 
in  other  models  the  residual  is  not  stationary. 

The  theorems  proved  above  jointly  with  the  theorems  in  works  [5,6] 
serve  as  substantiation  of  the  model  of  rhythmic  phenomenon  as  PKSP  and 
reveal  a  series  of  its  significant  properties,  and  also  show  that 
neither  additive  nor  multiplicative  models  can  be  considered  as  the 
most  general  for  rhythmic  phenomena,  biological  in  particular,  which 
both  by  isolation  of  mathematical  expectation  and  by  normalization  are 
not  reduced  to  stationary  processes. 


10 


REFERENCES 


1.  BapT.ien  .M.  C.  BseaeHHe  a  Teopxio  c.iyMafiHux  npoueccoa.  H/1,  M.,  1958. 

..  B  0  I'l  'I  H  Ill  II  H  K.  C..  il  p  a  r  a  H  51.  FI. —  B  kh.;  OrCop  h  neptaaia  HH^opMauiiH,  29. 
•HayKoaa  .ayMKa*,  K.,  1971. 

3.  Pe.ibipaHj  H.  .M.,  PaftKoa  H.  A.,  Ill  m  a  o  a  F.  E.  KoMM>^aTllBHllle  HopMiipoaaH- 
iiu't  Ko.ibua.  <>H3MaTni3.  M..  1960. 

j.  I  .laabiuiea  E.  F.— Teopiia  aepodTHOcrefi  h  ee  npMMeHemiH.  1963.  VIII.  2. 

^  i%9*  ^  ^  ^  **  nepeaaMa  HHi)>opMauiiii,  22.  fHayKoaa  AyMKa>,  K-. 

*j-  J  P^a  r  a  II  3,  n. —  B  kh.:  OtOop  h  nepeaawa  imifiopMauiiH,  30.  *HayKOBa  ay’MKa*,  K-, 

7.  p  a  r  a  H  3.  FI.,  .M  ii  x  a  A  ,1  o  a  c  k  ii  fl  B.  H..  3  a  o  p  c  k  ii  A  H.  H.—  B  kh.:  Ot6op  ii 
nepeaaia  iiHipopMauHti,  28.  «HayKOBa  ayMKa»,  K..  1971. 

6.  Koae.i  C.  .M.-  B  kh.:  TpyAhi  IV.  M..  1959. 

M  a  .1  a  X  o  a  H.  <t).iyKTyamiii  a  aBTOKO.ie6aTeabHbix  ciicieMax.  cHayna*.  M..  1968. 

10.  Flpoxopoa  K).  B..  P  o  3  a  H  o  B  KD.  A.  Teopim  BepoHTHOcieA.  cHayxa*.  M..  1967. 

1 1.  P  0  M  a  11  c  H  KO  <ti..  C  e  p  r  e  c  a  F.  A.  Bonpocu  npMK.iaAHOro  aHa.iHaa  c-iynaAiibix 
npoueccoa.  <Coa.  paatio*.  M..  1968. 

12.  C  o  .3  b  C  e  p  r  e  0  A. —  B  kh.:  BiioaoruMecKiie  lacw.  «.Miip>.  .M..  1964. 

•  1.  T  o  a  c  T  o  A  H..  E  .1  e  A  K.  C.  AxycTiiKa  oxeaMa.  «.Mhp».  M.,  1969. 

14.  Brier  G.  W. —  .Monthly  weather  Keview.  1965.  93.  2. 


11 


DISTKEEUnCW  UST 


DISTRIBUnCW  DIRECT  TO  RECIPIEOT 


ORGANIZATiaJ 


MICROFICHE 


BOSS  DIA/RIS-2FI  1 
C509  BALL0C509  BALLISTIC  RES  LAB  1 
C510  R&T  1ABS/A\'EAD0CM  1 
C513  ARRADCTM  1 
C535  AVRADOCWTSAROCM  1 
C539  TRASANA  1 
Q592  FSTC  4 
Q619  MSIC  REDSTCNE  1 
Q008  Nnc  1 
Q043  ATMIC-IS  1 
E051  HQ  USAF/INEr  1 
E404  AEDC/DOF  1 
E408  AfWL  1 
E410  ASEdC/IN  1 
E411  ASD/FTD/ITIA  1 
E429  SD/IND  1 
POOS  DOEyiSA/I»I  1 

poso  dA/oavAixysD  2 

lOSl  AFTT/US:  1 
P090  NSA/OB  1 
2206  PSL  1 


Microficiie  Nbr:  rrD93C00014S 
Fr[>-ID(RS)  T-0S62-92 


